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A kicking sequence of the atom optics kicked rotor at quantum resonance can be interpreted as a 
quantum random walk in momentum space. We show how to steer such a random walk by applying 
a random sequence of intensities and phases of the kicking lattice chosen according to a probability 
distribution. This distribution converts on average into the final momentum distribution of the 
kicked atoms. In particular, it is shown that a power-law distribution for the kicking strengths 
results in a Levy walk in momentum space and in a power-law with the same exponent in the 
averaged momentum distribution. Furthermore, we investigate the stability of our predictions in 
the context of a realistic experiment with Bose-Einstein condensates. 
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I. INTRODUCTION 

In many contexts of statistical descriptions in physics 
and mathematics, random walk models are very useful 
for qualitative and quantitative analysis Random 
walks are one of the simplest stochastic process and they 
represent the basic model of diffusion phenomena and 
non-deterministic motion. When the steps for the dif¬ 
fusing particle have finite variance and are uncorrelated, 
the evolutions obey Gaussian statistics and standard dif¬ 
fusion is observed. 

In transport and diffusion in complex systems, the ba¬ 
sic hypotheses underlying the laws of ordinary Brownian 
motion can often be violated. In particular, the lengths 
of the steps taken by the diffusing particles can have large 
fluctuations, and they can follow a probability distribu¬ 
tion with heavy power-law tails, obeying a generalized 
central limit theorem Hig. Examples of the so-called 
Levy a random-walk processes are observed in classic 
transport in complex materials [Mg and in many inter¬ 
disciplinary contexts in biology, ecology and economics 
niidn], making these processes a paradigm of trans¬ 
port and non deterministic evolution in the presence of 
large deviations. 

The concept of a classical random walk can be trans¬ 
lated into quantum random walks m using the entangle¬ 
ment between different degrees of freedom, e.g., between 
the spatial walk variable and a system intrinsic quantity 
such as spin |12j . An interesting question is how can a 
quantum walk turn into a classical one and vice versa. 
We address this question for a random walk in the mo¬ 
mentum space of kicked cold atoms. We will show how 
momentum distributions may be steered almost at will, 
realizing Gaussian diffusion and more complex a-stable 
distributions |3]. The idea of implementing random pro¬ 
cesses with scale-free behavior with cold atoms in optical 
lattices goes back to pioneering works from the quan¬ 
tum optics community, see, e.g., m which showed that. 


interestingly, classical random processes with power-law 
correlations can be mapped, under certain circumstances, 
onto a deterministic quantum problem (see also |14] for 
a more recent proposal in this context). 

In particular, it was proposed to exploit the quantum 
kicked rotor dynamics to realize directed quantum walks 
m and classical Levy walks by timing noise implying 
strong decoherence of the quantum motion |16j . As in 
this latter study, our proposal for a kicked atom realiza¬ 
tion of random walks misses the second degree of free¬ 
dom, which would be entangled with the momenta of the 
atoms. Yet, we will use the choice of a discrete phase of 
the kick potential to steer the walks, in close analogy to 
the coin degree of freedom in quantum walk theory m, 
which decides on the random direction of single steps of 
the walk. This is reminiscent of first implementations of 
quantum walks with classical optics where the coin was 
also no second degree of freedom but a random selection 
by optical elements such as beam splitters [I7|[IH]. In 
contrast to other studies mentioned above, our quantum 
motion is, in principle, fully reversible for a single re¬ 
alization of the walk, and classicality only results from 
the classical average over many realizations. Our results 
open the route to future investigations of our system, 
which fully include the second degree of freedom. This 
could be done for cold atoms and Bose-Einstein conden¬ 
sates using either internal states of the atoms (an effective 
spin nanniEo]) or more than one kick potential with dif¬ 
ferent phases and wavelengths to address independently 
different momentum classes m- 

Our paper is organized as follows: Sec. [IT] reviews how 
quantum resonant motion of the kicked rotor can be used 
to realize fast ballistic motion, either symmetric or di¬ 
rected in momentum space. Sec. |III| presents our central 
results showing how the input distribution of phases and, 
more importantly, of kick strengths convert into the ex¬ 
perimentally easily accessible averaged momentum dis¬ 
tributions at quantum resonance. Two specific cases are 
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discussed in detail: a Gaussian diffusive walk (III A I, 
which can be steered from the quantum to the classi- 
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cal regime, and heavy-tailed Levy walks (IIIB). Secs, 
and |V] underpin that the results of the previous section 
are indeed robust with respect to typical experimental 
limitations in atom-optics. We show this by checking the 
stability with respect to small detunings from resonance 
in sec. EYi and by taking into account additional compli¬ 
cations such as a finite window of feasible kick strengths 
and a finite width of a Bose-Einstein condensate in the 
Brilluoin zone of the periodic kick potential (i.e. in quasi¬ 
momentum) in sec. Sec. VI concludes the paper. 


II. QUANTUM WALKS AT QUANTUM 
RESONANCE 

Experiments on the quantum kicked rotor based on 
cold or ultracold atoms work with particles moving along 
a line periodically kicked in time by an optical lattice. 
Neglecting atom-atom interactions, the quantum dynam¬ 
ics are described by the following Hamiltonian in dimen¬ 
sionless variables (such that H = 1) 1^ : 

H{x,p,t) = — + kcos{a:)''^ S{t — Tt) . (1) 

tgz 

The kick period is r, the kick strength is A:, and T is a 
discrete time variable that counts the number of kicks. 
The periodic potential implies conservation of quasimo¬ 
mentum p. With the chosen units, P can take on allowed 
values between 0 and 1. Using Bloch theory, the atom 
dynamics from immediately before the {T — l)-th kick to 
immediately before the next T-th kick is then described 
by the one-cycle Eloquet operator [53]: 

Up,k{T) = e-i^(^+/5)V2 g-i/ccos(e) ^ (2) 

where M = —i^ is the (angular) momentum operator 
with periodic boundary conditions. In its usual realiza¬ 
tion with fixed kick strength, the full evolution over T 
kicks is thus described by 

Uj = Up,k{T)Up,k{T-l)...U0^k{2)U0,k{l). ( 3 ) 

A series of experimental investigations has looked at 
the so-called quantum resonant motion of the quantum 
kicked rotor [5H[5S] over the last decade, see e.g. and 
references therein. The interest in the resonant dynamics 
is mainly motivated by the type of ballistic motion with 
fast acceleration which can be obtained in this partic¬ 
ular parameter regime. We restrict our discussion here 
to the main quantum resonances occurring whenever the 
kick period t is not only commensurate to 27r but an 
integer multiple of it, i.e. r = 27r£, £ integer. Then, 
for specially chosen values of quasi-momentum, e.g. for 
P = 1/2 at £ = 1 and /? = 0,1/2 at A = 2, the energy of 
the /?-rotor asymptotically increases quadratically with 
the kick counter T [55] . From the theoretical point of 


view, in contrast to general values of the period, the main 
quantum resonances are accessible to analytical investi¬ 
gation [55H1II- Using the pseudoclassical approximation 
theory developed in [55H55] and reviewed in [53] , this re¬ 
mains true also for small detunings from resonant kick 
periods. 

Since the first factor standing for the free momentum 
evolution in ([^ is identical to one for the main quantum 
resonances at resonant quasimomentum, we can easily 
compute the momentum distribution obtained by just 
applying the second factor alone. The distribution after 
T kicks of strength k is therefore the same as the distri¬ 
bution after a single kick with strength kT. This gives a 
momentum distribution after T kicks for an initial state 
which is an eigenstate of momentum uq |26j 

P{n,T\no,k) = Jl_^^{kT). (4) 

Here Jm are ordinary Bessel functions with integer index 
m. The momentum distribution in the case of fc = 3 and 
no = 0 is plotted in Fig. [^a). We observe two domi¬ 
nant peaks, which move away from n = 0 at constant 
acceleration (i.e. we see ballistic motion). This linear 
increase in peak momentum is in accordance with the 
quadratic increase in energy occurring at quantum reso¬ 
nance [55H57] . Breaking the spatial-temporal symmetries 
of the kicked rotor, the quantum resonances permit the 
realization of directed transport as well (53] I5M33] . One 
example is shown by the thin dashed line in Fig. [^a), 
where the symmetry is broken by the initial state cho¬ 
sen to be a superposition of two momentum eigenstates 
|'!/'(n,0)) = (|n = 0)-|-e®‘^ \n = I))/-\/2 with relative phase 
(j) = 7r/2, see |531[3T| for details. 

We already see that the quantum resonances of the 
kicked rotor can be used to steer the evolution of the 
wave packet in momentum space. The symmetric ballis¬ 
tic peaks of Fig. 0^) are identical to the wave packets 
obtained from symmetric quantum random walks, as re¬ 
alized with optical pulses [niiig. Also, directed walks 
are possible when breaking the spatial-temporal symme¬ 
try by a correct choice of the initial state, just as in ran¬ 
dom quantum walks m- The random walk picture best 
applies in momentum space for our kicked system since 
each kick couples one specific momentum state to its 
neighbors, while during the free evolution the momenta 
are unchanged (“free flight”). The number of effectively 
coupled neighbors depends on the kick strength k. Prac¬ 
tically, we have just nearest neighbor coupling for fc ~ 1 
(see Eq. (|^ and |56|). In the following, we ask whether 
we can steer such a walk in momentum space by other 
means, e.g., by a random choice of kick strength and/or 
phases in the sinusoidal kicking potential. For this, we 
re-derive Eq. ([^ for these generalized cases in the next 
section. Introducing random couplings (whose strength 
changes with time) makes quantum evolutions classical 
due to decoherence, see e.g. refs. jUl [551 HSl I34H36] in 
the context of the quantum kicked rotor. Hence, we can 
also control the “quantumness” of the random walks in 
momentum space, not by changing our effective Planck 
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potential and/or the kick strength randomly from kick 
to kick. While for a fixed realization, i.e. fixing the se¬ 
ries {kt, 4’t}t=i,...,T, the evolution is fully coherent and 
quantum mechanical, averaging over many realizations, 
as above for various quasimomenta, makes the result in¬ 
coherent since the statistical average is a classical one. 
While this may seem to be a disadvantage, we shall see 
that this allows us to realize and implement in a realistic 
experiment with ultracold atoms, essentially any type of 
random walk we want. All that is required is to choose 
the appropriate distributions for our input parameters 
{kt, 

At quantum resonance with r = 2tt£, £ integer, and 
ignoring an n-independent phase, the Floquet operator 
simplifies to [50] 


FIG. 1. (color online), (a) Ballistic motion at quantum res¬ 
onance for a zero momentum plane wave initial state and 
k = 3, either symmetric around n = 0 (T = 30 with thin 
black solid line and T = 50 thick red/grey solid line), or di¬ 
rected at T = 70 for a superposition initial state, see text 
(dashed blue line), (b) Momentum distribution for a sym¬ 
metric walk as seen in (a) but a,t k — O.Stt, after T — 1000 
kicks averaged over lO'^ uniformly distributed quasimomenta 
P € [0, 27r). The red dashed line shows a fit according to the 
asymptotic formula AkKn^rP) derived in |26j . 


constant, but by classical averaging over different realiza¬ 
tions of the random walk. The effective Planck constants 
used here throughout are close to hes = r « 2Tr£ [25H27] , 
£ integer, and hence we are always in the deep quantum 
regime in this respect. 


III. LEVY WALKS IN MOMENTUM SPACE 


Up^kAT) = 


(5) 


with ^ = 7r£(2/3±l) mod(27r) and /3 the quasimomentum. 
The state after T kicks is given by applying consecutively 
the operator ([^ onto the initial state as shown in Eq. 
([^. Each of these applications is reversible. For exam¬ 
ple, a kick with kj and (j>j can be exactly reversed by 
a subsequent kick with fcj+i = kj and (j)j+i = 4>j + 
Hence, coherence is preserved during one complete kick¬ 
ing sequence. Such a time evolution leads to the following 
time-dependent wave function 


{e\'iljp,T) = (e 


— p-i^T cos( 6 I- 0 t) p-iiN -ifey-i cos{B-ij>T-l) 


_ _ g-*feicos(e-0i) 


Wrong choices of quasimomentum destroy the ballis¬ 
tic motion at quantum resonance. Typically, the larger 
the quasimomentum detuning from its resonant values 
the faster this will occur [26l • Hence, averaging quantum 
resonant motion over a distribution of quasimomentum, 
e.g., taken uniformly in the unit interval for simplicity, 
turns the averaged momentum distribution into an av¬ 
erage over peaks which stop at different stages in mo¬ 
mentum. The asymptotic result for T —>■ oo is then a 
(coarse-grained) momentum distribution in the integers 
n which decays like a power-law ^ This result was 
derived in |26| and rediscovered in m with some general¬ 
izations. Fig. j^b) shows realistic numerical simulations 
a,t k = O.Stt and T = 1000 which confirm the prediction. 
Therefore, averaging over various realizations of our ran¬ 
dom walk naturally leads to a momentum distribution 
with tails behaving like a Cauchy or Lorentzian distribu¬ 
tion. 

Here, we are interested in how to steer a random walk 
in momentum space by choosing the phases of the kick 


with 

T 

G{e, A {k,}, {</,}, T)=Y^ k, 008(0 - 4>, - (j - 1)0 
i=i 

= Re (e*® \Rt \= \Rt\ cos{0 + arg(i?T)) ( 7 ) 
and 

T 

Rt = Rt(A {k,}, {</,}) = E (8) 

i=i 

Fourier transforming ([^ into the momentum (n) repre- 
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sent at ion gives 


2ir 

{n\,Pp,T) = ^J e-™® (0|#,T) d0 
0 



27r 

J e-“® cos(e+arg(i?T)) ^^( 5 , _ 0) d9 

0 


gin arg(i?T) 

X 


27r 



—i{nO' -\-\Rt\ 


cos{0')) 


0 

-arg(i?T)-Te,O)d0', 


(9) 


where the substitution 9' = 9 + arg(i?r) is used. If the 
initial state of the atom is a plane wave with fixed mo¬ 
mentum pq =710+^0, then ^ takes the constant value 
^0 = 7'‘Z(2/3o ± l)mod(27r). For the plane wave initial 
state 





( 10 ) 


and neglecting all global phase factors that will cancel, 
the following probability distribution is obtained 


strength kes = \Rt\ (and (j)j = 0 for all j), or by a se¬ 
quence of Teff = I Rt I effective kicks with strength fc = 1 
{4>j = 0 for all j). This generalizes the discussion around 
Eq. which is regained for Rt — kT ai k = const. 
and ((> = 0. 

In short-hand notation we may now write 

OO 

P(n,T|no)= j d\RT\p{\RT\)Jl-nMT\)^ (14) 

0 


where the distribution p{\Rt\) represents the individ¬ 
ual averages over the distributions Pk{kj) and p,j>{4>j) as¬ 
sumed to be independent and identically distributed with 
given probability distributions for the kick strength and 
phase parts respectively. It is exactly the average from 
Eq. (141 in which we are mainly interested since it al¬ 


lows us to steer the random walk in momentum space by 
appropriate manipulation of /9(|i?T|)- 

Focusing on an initial state of the atom with momen¬ 
tum riQ = 0, the square of the Bessel function J^{x) in 


Eq. (14) can be roughly approximated for positive x and 
large pq as a Dirac delta function 


J^{x) « C5{x-\n\), 


(15) 


P(n,r|no,/3o,{fci},{<(>j}) = I(»t-I 

2 


h 


277 


-i{n-nQ)6 ^-i\RT\ cos{9) 




(11) 


where C is an appropriate proportionality constant. It is 
this property that is responsible for the ballistic “horns” 
seen in Fig. Ha). A numerical test and further justifica¬ 
tion for this approximation can be found in |39) . Within 
this approximation, the expression for the averaged mo¬ 
mentum distribution in Eq. (14) simplifies for large \n\ 
to 


In the second step, the definition of the Bessel function 
of first kind is used [5H] : 

Q+27r 

J„(z) = — / (12) 

27r J 

CK 

together with Jn{—z) = (—1)” J„(z) for integers n [55] . 

For simplicity, we first focus on the case of resonant 
quasimomenta, i.e. where ^ = 0, and refer to section 
|V 2| for the case of different quasimomenta. Then the 
argument function is 


P{n,T\no = 0) « / d\RT\ p{\Rt\) C S{\Rt\ - \n\) 

0 

= Cp{\RT\ = \n\). (16) 

Hence, the large scale behavior of the mean momentum 
is equivalent to that of the distribution p{\Rt\). Con¬ 
sequently, by choosing appropriate distributions Pk{kj) 
and p^{(j)j), p(|i?T|) and hence the large scale behavior 
of P(n,T|no = 0) can be steered in the desired way. 

A. Gaussian random walk 


T 

Rt = RrUkj}, {</.,}) = ^ k, , (13) 

1=1 

which corresponds to a random walk in the complex 
Argand plane for each realization of the parameters 
and {it'j}j=i,...,T- The length of each step 
is given by kj and the direction by the angles 4>j, |i?T| 
is the displacement from the origin after T steps. Hence, 
p{\Rt\) is the distribution of end displacements for the 
underlying walk. It is now easy to see that the momen¬ 
tum distribution after a kick sequence with parameters 
and is equivalent to a distribu¬ 

tion that would be obtained by a single kick with effective 


The simplest case of a random walk is obtained when 
only the phase of the kicking potential (j) is taken from 
a random distribution at fixed k. Since the phase is 
defined between zero and 27r there is not much choice 
for non-standard, e.g., heavy tailed distributions. We 
choose a Gaussian distribution, whose two limits are a 
delta function centered at some fixed value of the phase 
and a uniform distribution for a standard deviation much 
larger than 27r. The former case obviously reduces to the 
deterministic walk reviewed in Fig. Ha). This limit is 
highlighted in Fig. j^d). The latter case was analyzed in 
great detail in ref. |^[ , for random jumps in the quasimo¬ 
mentum. In a changed frame of reference this is identical 
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to jumps of the phase, the only difference lying in the 
fact that commutative shifts of quasimomentum induce 
additional shifts of momentum itself (an additional heat¬ 
ing effect not present here). We therefore can be brief in 
translating the phase walk in the complex plane into the 
language of a momentum walk. 

We may interpret the end position of the complex val¬ 
ued function Rt = iX,Y) as a vector in an Argand di¬ 
agram. The evolution of Rt corresponds to a diffusion 
in the complex plane, whose details are determined by 
the phase distribution. In all non-pathological cases, ac¬ 
cording to the central limit theorem, the distribution of 
(A, Y) for large T is hence given approximately by a two- 
dimensional normal distribution: 


f{RT = iX,Y)) 


1 1 r I 

-e ^ i i 


27T(7x o'y 


(17) 

where fj-i and (Ti {i = X,Y) &re the mean and standard 
deviation of X and Y, respectively. We obtain for the 
distribution of the random variable Rt that the mean 
E{Rt) = nx + i ^j|‘Y and the variance Var(i?T) = + 

a\. By considering the ^jS as taken from a symmetric 
distribution centered around zero, Yi{Rt) will be real and 
consequently /ry = 0. The smaller the width of such a 
distribution for the the more the random walk will 
be directed and E{Rt) = fix shifted along the positive 
real axis away from the origin. 


For a uniform distribution in [0,27r) for the it is: 





FIG. 2. (color online). Classical (a) to quantum (d) transi¬ 
tion for the random walk in the complex plane (for T = 50 
and k = 1). In (a) the (j>jS are distributed according to a 
uniform distribution in [0,27r), in (fe) — (d) the (j>jS are dis¬ 
tributed according to a normal distribution with mean zero 
and standard deviation a — 1.5, 1, and 0.1, respectively. In 
(d) we superimpose the deterministic result for a fixed phase 
(</) = 0) which is identical to the random phase result for small 
a, apart from deviations in the very center at small |n| < 5. 


E(i?r) = fj,x = 0 and = cr^ = T/c^/2, and we obtain 


KRt) = 


TrTkP' 


) — 


irTk'^ 


\Rt\^ 

Tk^ 


(18) 


Multiplying by 27r times the radius |i?T|, leads to the 
following distribution for the end displacement 


= ^exp 


\Rt\^ 

Tk^ 


\Rt\. (19) 


This distribution has a peak close to zero and shows a 
normal behavior. For non-uniform distributions of (pjs, 
o'x 7 ^ cry. Hence, Eq. 0 cannot be rewritten in a 
closed analytical form for the distribution p{\Rt\). The 
precise form of the phase distribution, p^{(j)j), can be 
used to steer the final distribution p{\Rt\), e.g. away 
from zero with px ^ 0. Because of Eq. (16), the 


averaged momentum distribution for large \n\ obeys 
P(n,T|no = 0) oc p{\Rt\), so that the final momentum 
distribution follows the distribution of p(|i?T|)- As seen 
in Fig. we may steer the random walk from quan¬ 
tum, panel (d), to classical, panel (a), by controlling the 
widths ax = ay of the phase distribution. 


B. Power-law walks in momentum space 


We have seen above that random phase shifts alone 
cannot lead to more interesting non-Gaussian distribu¬ 
tions. In consequence, we must include the possibility of 
randomizing the kicking strength k in Eq. (13) as well. 
Obviously, a Gaussian distribution with fixed mean and 
standard deviation would have the same consequences as 
the case just discussed in the previous subsection (with 
the only difference being a change in the final Gaussian 
momentum distribution). 

However, more interesting choices are possible since 
k is not bounded from above, just from below at fc = 
0. For simplicity, we can interchange a distribution of 
just positive values of k with a symmetric distributions 
around zero with. This is possible since the phase (jjj may 
be chosen from a discrete uniform distribution, which 
only takes two values with equal probability p = 1/2: 


= Y ^ !■ (20) 

P=2 

From the form of the kick evolution operator, Eq. (§, it 
can be seen that a tt phase shift is equivalent to a sign 
change in front of k. Such a phase shift is experimen¬ 
tally feasible by varying the relative position between 
the atoms and the kicking standing wave [32l I4nil42j . 
Furthermore, asking which of the two phases should be 
chosen before each kick, is analogous to the coin toss in 
standard quantum random walk algorithms, see e.g. [12j . 

In the following, we focus on so called a-stable distri¬ 
butions S'(a, 7 ,/r) for the probability distribution of the 
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kick strength. These have several advantages, the most 
important one is that they are closed under convolution 
for a fixed value of the parameter a, in the sense that 
the sum of N independent and identically distributed 
random variables Xj ^ 5 '(a, 7 ,/r) is again distributed 
with ^ S{a, Nj, Nfj,) [3]. Here the first pa¬ 

rameter a G ( 0 , 2 ] is called the characteristic exponent, 
it describes the tail of the distribution. 7 > 0 is a scale 
parameter characterizing the width of the central part of 
the distribution and /x G K determines the position of its 
center [3]. For a < 2, the stable distributions have an in¬ 
finite variance and asymptotically decrease as 
[3]. Hence, by choosing /x = 0, the resulting distribution 
S'(a, 7 , 0 ) is a distribution symmetric around zero with 
power-law tails. 

Most of our data is produced now when the kjS obey 
a Cauchy distribution, as a special case with a = 1. We 
formally assume that k can be negative to simplify the 
argument. Here an explicit form of = 0) can be 

given: 

Pk{kj) = - 2 ,2 ■ (21) 

Let us also assume for a moment that (j) = const, (in¬ 
cluding the case of jumping between 0 and tt for the real¬ 
ization of effectively negative k). From the fact that the 
kick strengths kj are Cauchy distributed with 5 '( 1 , 7 , 0 ), 
it follows that the sum Rt = exp(—i(/)) kj is also 

Cauchy distributed with Rt ^ S'( 1 ,T 7 ,0). |i? 7 ’| is then 
the “length” of the end positions of a one-dimensional 
Levy walk on the real axis. The distribution of end dis¬ 
placements is therefore given by 

PURtD = 



FIG. 3. (color online). Numerically obtained averaged mo¬ 
mentum distributions with 50,000 realizations for the biS and 
kjS at T = 50, no = 0 and ^ = 0. For the kjS a Cauchy dis¬ 
tribution with scale parameter 7 = 0.5 is used taking into 
account just positive values of k (see also Eg. ( |30[ ) below), for 
the (f)jS the discrete distribution from Eg. (20 1 (blue line with 
dots) and a uniform distribution (red dotted line with aster¬ 
isks) is considered. The thick black line represents a power 
law with exponent 2 (or an a parameter a = 1 ) fitted to the 
data. 


results are identical in the tails, showing a power law with 
exponent a = 1 at large momenta |n| > 30. 

To summarize by choosing the kicking strengths from 
an appropriate power-law distribution pk{kj), a Levy 
walk can be realized in the momentum space of the atoms 
with the same asymptotic power law as the one of the in¬ 
put distribution pkikj). The momentum distribution is 
the quantity typically measured in an atom-optics exper¬ 
iment [221 [23] . 


IV. STABILITY WITH RESPECT TO 
DETUNINGS: e-CLASSICAL ANALYSIS 


For the large time asymptotics, it does not matter 
whether (p = const, as just assumed (corresponding to 
an effectively one-dimensional walk) or drawn from a uni¬ 
form distribution in [ 0 , 27r) (fully two-dimensional walk in 
the complex plane), or anything intermediate. For large 
|i?T|, the distribution will always be similar to 5'(1,7,0), 
with a scale factor 7 < T 7 since the optimal spread is 
achieved exactly for the one-dimensional walk. Thus the 
tails are power-law distributed in all cases, with the same 
power-law (i.e. the same a) as the kick strength distri¬ 
bution. 

Figure]^ shows in a log-log plot the results of numerical 
simulations carried out at quantum resonance for 50,000 
realizations of the walk for with different series and 

{fcj}s at T = 50, no = 0 and ^ = 0 (resonant quasimo¬ 
mentum). The plot shows a comparison of the averaged 
momentum distribution obtained by considering for the 
kjS a Cauchy distribution with only positive values, scale 
parameter 7 = 0.5, and for the (^js either 0 or tt (blue 
line with dots) or a uniform distribution in [0, 27r) (red 
line with asterisks). Recall the possibility of (j)j > tt ef¬ 
fectively corresponds to negative kicking strength. Both 


Quantum mechanically, there is no general analytic so¬ 
lution outside the quantum resonance condition for the 
quantum kicked rotor dynamics, and we may refer only 
to numerical simulations, as we will do below. However, 
if we detune the kick period r just slightly from the res¬ 
onance conditions, we may still be able to estimate the 
true quantum motion using the pseudoclassical model in¬ 
troduced in [2SI2II- In this model, the absolute value of 
the detuning |e| plays the role of Planck’s constant, and 
hence the theory has a semiclassical limit at exact quan¬ 
tum resonance. The e-classical dynamics are described 
by the following discrete map which relates the variables 
I and 9 from the j-th kick to the next one: 

Ij+i = Ij + kj^i sin(0j -|- (/fj+i) /23'\ 

= 6j + sgn(e)/j_|_i mod 27r ' ^ ^ 

Here, I = J-|-sgn(e)[ 7 r£-|-T/ 3 ] and J = \c\p is the rescaled 
momentum. Similarly, a new effective kick strength 
is defined by kj = \e\kj^ which is multiplied by the 
small detuning in this model, making the classical phase 
space nearly integrable. Provided that the effective kick 
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strength k remains small, the map (231 yields a good ap¬ 
proximation of the true quantum motion [261 ET]. The 
momentum distribution studied in the previous sections 
corresponds to a distribution of p{T,9o) = J{lT,So)/\e\ 
for the initial ensemble of 9o (uniformly distributed in 
[0,27r) for modeling a plane wave initial condition with 
fixed momentum) and the realizations of the random 
variables {kj}. The momenta It after T kicks are 
given by 


T-l 

It = lo+ '^kj+ism{9j + (l>j+i). (24) 

j=o 



To begin with, we look at the case of fixed k and broadly 
distributed phases, just as in subsection |mg Then the 
single momentum changes are uncorrelated and the mo¬ 
tion will be diffusive just as assumed for the determin¬ 
istic standard map in the chaotic regime |43l |44| . For¬ 
mally, the momentum increase AIt = It — Iq is related 
to our distribution function Rt from Eq. (13), by using 
= ~(^j + 4‘j+i) s-nd rewriting 


AIt = — Im 



f'j 1 — _ 


T 

i=i 




(25) 

The random walk obtained in the variable I is then just 
a projection of the walk in the complex plane onto the 
imaginary axis. Retranslating into the original momen¬ 
tum variable p, gives a result very similar to the case 
studied in sections III A and [III B[ when the phases (j)i 
strongly vary in the interval [0, 27r). 


For a fixed phase (p, the momentum increase will dom¬ 
inantly grow by the sum of kick strengths, since we may 
assume a random reshuffling of the angles due to the ran¬ 
dom change in k. This is similar to the argument used 
to derive momentum diffusion in the deterministic stan¬ 
dard map |43l 04] , where chaos there is substituted here 
by true randomness. In other word, we may approximate 


T _ T 

kj sin(6»j +^) - var (sin(0^- -f <)>)) ^ kj 
i=i i=i 




i=i 


(26) 

with var {siii{9j + 4>)) « 1/2 from the approxima¬ 

tions (sin(0j+i) sin(0j)) « 0 and (sin(0j) sin(0j)) « 
Y dxsin^x = 1/2. Hence, we expect that the walk 
in momentum space can be steered by the distribution 
of kick strengths close to quantum resonance making it 
robust with respect to experimental imperfections. Our 
expectation, based on the crude approximation from Eq. 
( |26| ), is confirmed by full quantum simulations at finite 
detunings, which are shown in Fig. [^ as well as by 
pseudo-classical simulations based on the above reason¬ 
ing (not shown here). 


FIG. 4. (color online). Averaged momentnm distributions 
computed by evolving the quantnm map generalizing Eq. (§ 
with randomly distributed kj from a Cauchy distribution (no 
phase changes). Parameters are r = drr -|- e, 7 = 0.5 after 
T = 50 kicks, and 20000 realizations are taken into account. 
Data are shown for three different detunings e = 0.001 (black 
dotted line with circles), 0.01 (red/grey thick solid line), and 
0.1 (blue dashed line). The different curves only differ in the 
center close to zero momentum while the tails are identical. 
The case e = 0.001 is not distinguishable from the result at 
exact quantum resonance, see blue line with circles in Fig. 
The black solid line shows a power-law ht with exponent 2 as 
expected for the Cauchy case. 


V. EXPERIMENTAL IMPLEMENTATION 

We now investigate how the predictions of the pre¬ 
vious section could be applied in a state-of-the-art ex¬ 
periment with Bose-Einstein condensates, such as the 
one successfully running at Oklahoma State University 
[31 iS ED 0107]. We focus, in particular, on pos¬ 
sible problems in the implementation of random walks 
with heavy tail statistics. It is clear that the experimen¬ 
tal challenge is two-fold: first, relatively large momenta 
should be observed, at least over one or two orders of 
magnitude. Second, the signal in the wings should be 
large enough to give an experimental signal above the 
noise level, see e.g. 0105]. Next, practical problems in 
the realization of the ideal walk must be considered, e.g., 
the challenge of rapidly changing the kick strength from 
kick to kick over an appreciable range of values, or the 
finite width of a Bose-Einstein condensate in quasimo¬ 
mentum. In the subsequent subsections we will address 
these and other possible effects. 


1. Fixing the kicking strength for a seguence 

The first issue to address is that it may be easier not 
to change the kick strength for every kick during a mea¬ 
surement sequence. Therefore, our proposed idea is to 
fix the kick strength kj = k for j = l,...,r during a 
sequence, and to choose for every repetition of the ex- 











periment a different k drawn, for example, from a heavy 
tailed power-law distribution. Finally, an average over all 
results can be performed to obtain the walk momentum 
distribution. In such a case, the individual experimen¬ 
tal run is no longer a random walk in momentum space 
since k is fixed. However, if the potential shifts (j)j can 
be adjusted precisely during a sequence, a Gaussian walk 
trajectory, as studied in |III A[ can be realized for each 
run. We will show now that in this case, the final aver¬ 
aging over many realizations eventually leads to a heavy 
tailed momentum distribution as well even if k is fixed 
for each repetition but chosen randomly from an a-stable 
distribution. 


Following Eq. (141, at quantum resonance, the mean 
momentum distribution after T kicks is now explicitly 
given by 


P{n,T\no) = J dk pkik) J dcjyi p^{(l)i) ... (27) 


d(j>T Pp{4>t) Jn-no (I^t|) 


where only one integral over k needs to be considered. 
We must now derive an expression for the distribution 
p{\Rt\), with Rt = ke~^’^p a random walk in the 

complex plane with fixed jump length k for each indi¬ 
vidual repetition of the experiment. We focus on the 
simplest case in which the (j)jS are chosen from a uniform 
distribution in [0, 27r). From section III A it is known 
that for fixed k and uniformly distributed d<jS, the distri¬ 
bution of end displacements is given by Eq. (19), which 
leads to the following mean displacement 


\RT{k)\ = I\RTik)\p{\RT{k)\)d\RT{k)\ 
0 



(28) 

The inflnence of the c^^s, which determine the displace¬ 
ment of a single sequence with fixed k, is now neglected 
since it is small compared to the possibly large jumps in 
k. Then we can assume that the displacement of a single 
sequence is given approximately by the mean 


\R^{k)\^\RT{k)\ = ^\k\. (29) 

Since this mean displacement is linearly dependent on 
|fc|, the resulting distribution p{\Rt\) is again completely 
dominated by the chosen distribution for |fc|. To be spe¬ 
cific we now consider a “positive” Cauchy distribution 
with fc > 0 


p{k) 


27 

7r(72 -I- fc2) ’ 


(30) 


where the factor two corrects the normalization due to 
taking into account only positiv values. By substitution 




FIG. 5. (color online) Log-log plots of the numerically ob¬ 
tained averaged momentum distributions (blue symbols) via 
simulations with 50,000 realizations for the and k,T = 50, 
no = 0 and C — 0 - For the cpjS a uniform distribution is con¬ 
sidered, for k the Cauchy distribution Eq. ( |30| l with scale 
parameter 7 = 0.5 [panel (a)] and 7 = 0.1 [panel (b)]. Ad¬ 
ditionally, we show the averaged momentum distribution ob¬ 
tained by computing Eq. ( |14| | with the derived analytical 
expressions for p{\Rt\) (black solid line) from Eq. (|31[). 


p(fc) dk « p{\Rt\) dlRxl we arrive at 


p{\Rt\) 





TT I 7j + I^Tp 


(31) 


As in section IIIB the distribution of fc’s from Eq. (30) 


carries over to a Cauchy distribution for the end dis¬ 
placement |i?T|, although with a smaller width, 
as compared to T 7 . Eigure confirms the validity of 
this approximate result. The plots also show that, at 
smaller 7 , the power law starts at lower values of mo¬ 
mentum, while at larger 7 the power-law starts at larger 
momenta. 


2. Taking into account the quasimomenta 

Up to now, always resonant quasimomenta have been 
considered, corresponding to ^ = 0 in Eq. (|^. Modern 
BEC experiments allow one a relatively precise adjust¬ 
ment of the /3 distribution which is approximated to be 
Gaussian [IH |46] . Sticking to the quantum resonance at 
T = 27r (i.e. f = 1), we consider the distribution centered 
around the resonance value /3 = 0.5: 

1 {( 3 - 0 . 5)2 

Pp{P) = e 2 .:^ , (32) 

(7V27r 

with a = i//2/-\/ln(2) 2, where v k. 0.05 is the full width 
at half maximum (FWHM) measured in [JT]. An ade¬ 
quate value for the number of considered /3s is H = 10,000 
corresponding to the typical atom number in the exper¬ 
iment [33l l4nil4^ l45tiT7] . This implies the assumption 
that each atom has its own quasimomentum, neglecting 
interactions, and hence we assume that the initial atomic 
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ensemble is an incoherent mixture of plane waves with 
different but fixed quasimomenta. 

Hence, given the probability distribution after a 
kicking sequence for an atom with fixed (3 = /3o, 
P(n, T|no,/3o, {fcj}, the “real” momentum distri¬ 

bution, under the incoherent approximation, is given by 
the following classical average 


! dppp{j3) P{n,T\nQ,l3,{k^},{(j)j}). (33) 

0 


At 

P(n, 

( 0 , 


quantum resonance and for a 

T\no,P,{kj},{(j)j}) = 

where Rt now takes the extended form 


hxed 

see 


P, 

Eq. 



Rt = Pt({A:,}, {h},m) = E ( 34 ) 


i=i 


This formula underlines the similarity between ^ (and 
hence p) and pj acting both as effective potential dis¬ 
placements. Both have an influence on the direction of 
the single steps performed in the underlying walk in the 
complex plane. If we consider a case in which the pjS 
are uniformly distributed in [0,27r), the potential is com¬ 
pletely randomly shifted in 9 space for every kick. Fur¬ 
ther shifts arising from the different P are therefore av¬ 
eraged out and negligible. Thus, in this special case, the 
effect of the /3s can be absorbed in the phase average. 

Figure shows the results of numerical simulations 
for which the additional averaging over the /3s is done. 
We consider the case in which the pjS are uniformly dis¬ 
tributed and k is fixed and chosen from the Cauchy distri¬ 
bution (301 for each repetition. First, for the simulation 
in panel^), 50,000 realizations for the pjs (j = 1, 
and k are considered at only one non-resonant value for 
/3, i.e. at T = 27r with /3 7 ^ 0.5. In all cases, one obtains 
a similar result and even for P = 0, which is the value 
the furthest from resonance, no significant difference is 
observed. In panel ( 6 ), several normal distributions for 
P are considered with a smaller number of realizations. 
Again, for all considered values of iz, no signihcant differ¬ 
ence is seen. The differences for higher |n| result from the 
fact that even R = 25,000 realizations are not sufficient 
to obtain always a power law up to the tails at |n| > 10 ^ 
with signihcant statistics. In the next sections, simula¬ 
tions with smaller values for R are done, thus giving the 
opportunity for additional checks with a more realistic 
number of quasimomenta. 



FIG. 6. Log-log plots of numerically obtained averaged mo¬ 
mentum distributions via simulations which take into ac¬ 
count the quasimomenta. We consider the </)jS uniformly dis¬ 
tributed, k fixed for each run (see section V 11 and chosen 
from the Cauchy distribution ( |30[ ) for each repetition; the pa¬ 
rameters are T = 50, no = 0 and r = 2^. In (a) the results 
are plotted by considering R = 50,000 realizations for the pjS 
and fc and the different single values for /3 given in the leg¬ 
end. In (b) the results are plotted by considering R — 25,000 
and B = 200 values for P chosen from a normal distribution 
centered around the resonance value /3 = 0.5 with different 
values for the FWHM u, see legend. 


laboratory, one may realize without problems a range 
of k between about 0.1 and 5. The lower limit is set by 
experimental uncertainties, which make it hard to distin¬ 
guish between k = 0 and fc = 0.1 kicking signals. Hence, 
only a small part of the Cauchy distribution for k > 0 
can be considered: 


3. Limits for the kicking strength and the number of 
repetitions 


Pk{k) = 


const. I 

0 


for A:niin P: k pi A:niax 
otherwise 


(35) 


Another experimental problem is that no infinite range 
for the kicking strength can be realized, which would be 
necessary to cover the entire tail of a power-law input 
distribution. Given the typical laser intensities in the 


where const, is an appropriate normalization constant. It 
is suitable to use 7 = 0.1 since for such a small scale pa¬ 
rameter, the power law sets in early and covers the acces¬ 
sible interval for k. Furthermore, since one measurement 
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takes about 30 seconds, 1,000 repetitions are reasonable. 
Because of the expected cutoff at large momenta, it is 
also sufficient to consider only 1,000 realizations from a 
statistical perspective. 

Figure shows the results of numerical simulations 
for which all these realistic experimental parameters are 
taken into account. For (3, B = 10,000 realizations and 
a normal distribution with v = 0.06 is considered. As 
above, for the (j)jS a uniform distribution is assumed and 
k is chosen from a Cauchy distribution with the above 
limits for each repetition. We consider a small upper 
limit /cmax = 5 with fcmi„ = 0 in panel (a). One observes 
that the limit for the kicking strength leads to a cut-off 
in the averaged momentum distribution. This is easily 
explained with the delta approximation for the square of 
the Bessel function in Eq. (15). Since k is limited, the 
values for |i? 7 ’| are limited and hence the momentum dis¬ 
tribution is nearly zero for all |n| > |i?T|max- The cut-off 
value in the momentum distribution depends linearly on 
the maximum limit for k, this is in accordance with the 
mean displacement, which is also linearly dependent on 
k. We verified numerically that there is indeed a lin¬ 
ear correlation between the cut-off value in the momen¬ 
tum distribution and the maximum kicking strength (not 
shown here, see [39] for details). In panel ( 6 ), the result 
of exactly the same simulation as done for (a) is shown, 
with the only difference that here also the minimum limit 
fcmin = 0.1 is taken into account. There are small differ¬ 
ences between the resulting distributions. In the second 
case, the power law sets in slightly later. Furthermore, 
in the second case the power law region is slightly shifted 
upwards in the vertical direction. The reason for this ef¬ 
fect is that by considering a minimum limit, an interval 
of the Cauchy distribution (from 0 up to 0.1) is cut which 
has a relatively large probability. Hence, due to the new 
normalization, the values for k in the range from 0.1 up 
to 5 have a higher probability compared to the first case. 
In other words, the constant in front of the distribution 
in Eq. (35) depends on the specific limits. Consequently, 
the momentum values resulting from these values of k 
are more likely. 


4- Fluctuations in the kicking strength and optimizing the 
signal-to-noise ratio 


Another experimental problem is that k is not com¬ 
pletely constant during a sequence because of two inde¬ 
pendent experimental artifacts. The first one is that not 
all the atoms sit at the same spot and hence they may 
feel a slightly different intensity of the kick potential. 
Additionally, there are time-dependent drifts in the ex¬ 
periment, especially when data is taken over a relatively 
long period. Considering such fluctuations are of the or¬ 
der Ak = ±0.1, consistent with the lower cutoff chosen in 
the previous subsection, one does not observe significant 
differences in the results, see ref. [331 for details. 

In previous atom-optics experiments weak spontaneous 



FIG. 7. Log-log plots of the numerically obtained aver¬ 
aged momentum distributions (blue symbols) via simulations 
which are close to a real experiment. For /3, 10,000 values 
taken from a normal distribution with FWHM iz = 0.06 are 
considered. For the <j)jS a uniform distribution is regarded 
and k is chosen from a limited Cauchy distribution for each 
repetition. The used limits are femax = 5 and in (a) fcmin = 0 
and in (6) fcmin = 0.1. R = 1,000 realizations are considered 
at T = 50 for no = 0. Both panels also show a power law with 
exponent 2, corresponding to a Cauchy distribution, fitted to 
the data (black solid line). 


emission was changing the quasimomentum of the atoms 
and also heating them slightly (231 . While this is 

a nuisance, for instance, for the observation of dynamical 
tunneling it should not be a problem in the present 
context. Changing randomly quasimomentum effectively 
leads to a broader distribution of quasimomenta, and it is 
similar to changing randomly the phase of the kick poten¬ 
tial (see discussion in sec. Ill A). Since both of these ef¬ 
fects do not alter very much our predictions, spontaneous 
emission should not hinder the observation of power-law 
tailed momentum distributions. 


The final problem we consider here might be noise in 
the detection procedure (e.g. due to pixel fluctuation in 
the CCD camera at low contrast), which makes it diffi¬ 
cult to see small signals in the momentum distribution for 
large |n|. Since due to all previous considered limitations, 
the power law obtained in the momentum distribution 
only covers a small momentum range and for this range 
the probability ratio P{n = 0)niax to P{n)^in is about 
three orders of magnitude, this might not be a problem. 
However, a possible way to avoid the noise problem is 
to choose for the k distribution a power law that de¬ 
cays slower than a Cauchy distribution, i.e. 0 < a < 1. 
Hence, the ratio F’(n)max to P(n)n,i„ will not grow too 
fast, allowing for a better experimental signal. Together 
with an appropriate choice of the width parameter 7 , see 
discussion at the end of subsection |V ![ this should give 
the possibility for the experimental implementation of 
steered random walks with a Bose-Einstein condensate. 
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VI. CONCLUSIONS 


We have shown how it is possible to steer the aver¬ 
age momentum distribution of kicked ultracold atoms by 
applying random sequences of kicks with specific distri¬ 
butions of phases and kick strengths. 


Our central result is contained in formulae (14) and 


(16), respectively. They show that the chosen distribu¬ 


tion of kick strengths converts into a momentum distri¬ 
bution with the same asymptotic scaling. This allows 
for future experimental investigations of complex random 
walks, e.g., power-law distributed Levy walks in momen¬ 
tum space. Random walks are an important part of many 
biological, social, and physical systems [D12 H E HU]. 
Having a robust scheme for the implementation of com¬ 
plex classical walks with great stability will enable new 
insights into many non-deterministic transport processes 
in nature with large fluctuations. The robustness of 
our prediction with respect to experimental limitations is 
confirmed by our detailed analysis of the effects of small 
detunings from quantum resonance, of limits in the kick- 
strength distribution, and of a finite quasimomentum dis¬ 
tribution of the atoms (see secs. IV and 0 . 

A random choice of the phase of the kicking lattice 


from the two values 0 and tt mimics the coin toss of a 
quantum walk since it effectively changes the sign of the 
kick potential, and hence determines the direction of a 
single kick or step of the walk. Quantum random walks 
have been realized in different setups with single parti¬ 
cle control [50H53] . Extending our study to include the 
necessary additional degree of freedom - to control the 
single step of the walk simultaneously entangled with the 
center-of-mass momentum degree of freedom of the atoms 
- may permit one the implementation of stable quantum 
walks with Bose-Einstein condensates with much more 
particles. Moreover, in such improved walks with ultra¬ 
cold atoms, the quantum to classical transition [Ml IM] 
in the walk behavior may be investigated as partly an¬ 
ticipated by our Fig. 
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